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Challenges:  Reactive  LES 
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Algorithm  comparisons: 

-  identical  subgrid  modeling 

-  differences  reside  in  numerics 
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Ref:  2013  -  Cocks,  Sankaran,  Soteriou,  “Is  LES  of  reacting  flow  predictive?  Part  1:  Impact  of  Numerics” 


Need  to  determine  BEST  discretization  schemes  for  Reacting  LES 
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Objectives 


GOAL: 

damp  high  frequency  errors  while  preserving  low  wave 
content  (ie:  low-pass  response) 

1)  compare  damping  character  of  Artificial  Dissipation  and  Filtering 

2)  formulate  filter  as  an  equivalent  Artificial  Dissipation  scheme 

-  consequence  of  filter  damping  for  stiff  problems 

3)  insight  on  achieving  “ideal”  low-pass  response  for  general 
problems 


•  von  Neumann  Analysis 

•  Crank-Nicolson  w/  6th  order  central  differencing 
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von  Neumann  Analysis 

1 D  Euler  System  (quasi-linear  form): 

dQ  a  dQ  „ 

— +  A— =  0 

dt  dx 


Q,=^Q{k)eib: 


QM  =  Y,Q(k)  gik(x+Ax)  _  Y,Q(k)  gikAx^ikx  kAx  =  0 

k  k  with  0  €=  [-;r,;r] 


Q"+'  =  G(e)Q 


n 


Eigenvalues  of  the  amplification  matrix  specify 
growth  factor  and  phase  errors. 
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Artificial  Dissipation 


dQ 

dt 


^(-ir  (AxT^e^ 


d2mQ 

dx2'" 


governing  equations  augmented  with  dissipation  terms 


e2  =  1  /  2 
e4  =1/12 
e6  =  1  /  60 


upwind  biased  stencils 


Filtering  (Explicit) 


2m 


\2m 


m 


dx 


2m 


smoothing  of  solution  as  post-process  of  integration  step 


represent  total  amplification  of  filter  scheme  as: 


G(d)  =  R(d)G*(d) 


with 


|«(9J|2<i 
/?(0  =  ±jt)  =  O 


NOTE:  filter  is  purely 
dissipative  and  does  not 
alter  original  scheme’s 
phase  behavior 
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Filtering  (Implicit) 


2m 


\2m 


m 


dx 


2m 


Q,= 


l  +  X  S2m(&xf- 


1 2  m 


m 


dx 


2m 


smoothing  of  solution  as  post-process  of  integration  step 


represent  total  amplification  of  filter  scheme  as: 


G(d)  =  R(d)G*(d) 


with 


|«(9J|2<i 
/?(0  =  ±jt)  =  O 


NOTE:  filter  is  purely 
dissipative  and  does  not 
alter  original  scheme’s 
phase  behavior 
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Artificial  Dissipation 
Growth  Factor 
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Explicit  Filter:  Growth  Factor 


Shapiro  Filter  (1975) 


Q,= 


l  +  Sm(Axfm 


\2m 


dx 


2m 


Q*  with  5 


8 


Distribution  A- Approved  for  public  release;  distribution  is  unlimited 


2nd  Order  Implicit  Filter: 
Growth  Factor 
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2nd  vs.  4th  Order  Implicit  Filters: 

Growth  Factor  Error 


4th  order  Lele 
Filter  (1992) 


[1+5(A*)2i] 
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2nd  Order  impFil  (del  =  1e-1) 
2nd  Order  impFil  (del  =  1e-2) 
2nd  Order  impFil  (del  =  1e-3) 
2nd  Order  impFil  (del  =  1e-4) 
2nd  Order  impFil  (del  =  1e-5) 
4th  Order  impFil  (del  =  1e-1) 
4th  Order  impFil  (del  =  1e-2) 
4th  Order  impFil  (del  =  1e-3) 
4th  Order  impFil  (del  =  1e-4) 
4th  Order  impFil  (del  =  1e-5) 


—  -  -  500th  Order  expFil 
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Filtering  as  an  Artificial 
Dissipation  Scheme 
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dispersive  terms: 

•  restore  phase  of  original  scheme 
dissipation  term  scales  as  e2  ~  1  / CFLu+c 

•  increased  damping  w/  decreasing  time-step 
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Effect  of  CFL 


*  AD  shown  is  6th  order 


AD  (CFL  =  I  e-2) 
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Cumulative  Low  CFL  Damping 


CFLu+c  =  10-4 
Nsteps  =  104 
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Summary 


•  Filtering  is  a  form  of  artificial  dissipation 

-  damping  behavior  more  predictable  and  tunable 

•  Explicit  filters  require  very  high  order  for  low-pass  response 

-  overly  dissipative  for  small  time-steps 

•  Implicit  filters  can  be  efficiently  designed  for  low-pass  response 

-  superior  to  artificial  dissipation  or  explicit  filters 
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Staggered  Grid 
Von  Neumann  Analysis 


d  Q  dE 
dt  dx 


=  0 


ID  Euler  Eqns 


Eigenvalues  of  the  amplification  matrix  specify 
growth  factor  and  phase  errors. 


Qn+1  =  GQ 
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Staggered  Grid  Scheme/  Quasi-Linear  Form 
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Runge  Kutta  4:  von  Neumann 


Collocated  Staggered 


Kinetic  Energy  Preservation  (KEP) 


“in  computations  of  turbulent  flow  fields,  dissipative  errors  show  up  at  the  level  of  kinetic  energy”  (Mahesh  2004) 
Robust  at  inviscid  limit  (Re  00) 


Incompressible  Flow: 
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Compressible  Flow: 


K  =  14  U|2  bounded  and  constant  at  inviscid  limit 
KEP  schemes  satisfy  secondary  equation  discretely 
Richtmeyer  &  Morton  (1967) 

Arakawa  (1966) 
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Discrete  analogue  seeks: 

Accurate  transport  of  KE  accurate  physical  transfer  of  energy:  E  =  KE  +  Uint 
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KEP:  Applied  to  ID  Euler 


(Collocated  Grid)  | 


compare  Crank-Nicolson  (CN)  with  Fully  KEP  scheme  (F-KEP) 
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discrete  secondary  equation  satisfied  to  machine  zero  if  KEP 
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Evaluating  KEP:  u 


Mach  =  0.859 


disturb  =  0.01% 


-  F-KEP  always  secondary  conservative 

-  CN  is  KEP  for  low  compressibility  effects 
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j.  n-7  Crank-Nicolson  Scheme 
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